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Abstract 

We discuss the thermodynamics of the Af = 2*, SU(N) gauge theory at large 't Hooft 
coupling. The tool we use is the non-extremal deformation of the supergravity so- 
lution of Pilch and Warner (PW) [hep-th/0004063], dual to Af = 4, SU(N) gauge 
theory softly broken to Af = 2. We construct the exact non-extremal solution in five- 
dimensional gauged supergravity and further uplift it to ten dimensions. Turning to 
the thermodynamics, we analytically compute the leading correction in m/T to the free 
energy of the non-extremal D3 branes due to the PW mass deformation, and find that 
it is positive. We also demonstrate that the mass deformation of the non-extremal D3 
brane geometry induces a temperature dependent gaugino condensate. We find that 
the standard procedure of extracting the Af = 2* gauge theory thermodynamic quanti- 
ties from the dual supergravity leads to a violation of the first law of thermodynamics. 
We speculate on a possible resolution of this paradox. 
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1 Introduction 



Over the last few years, gauge theory/string theory duality [1] (see [2] for a review) 
has proven to be a very useful tool to address nonperturbative questions in gauge 
theories. Essentially, this duality states that, say, four dimensional gauge theories 
at large 't Hooft coupling, gy M N 3> 1, can be described by dual string theories in 
weakly curved supergravity backgrounds. In the large N but fixed 't Hooft coupling 
limit, the string coupling vanishes, and one can consistently restrict the string theory 
side of the correspondence to the massless sector of type IIB supergravity. Thus the 
computations on the supergravity side shed light on the nonperturbative gauge theory 
dynamics. On the other hand, nonperturbative effects in the gauge theory potentially 
could tell us something new about the dual supergravity. This is indeed the case 
with finite temperature phase transitions in gauge theories [3-7]. Specifically, in [3], 
it was demonstrated how the (kinematic) confinement-deconfinement phase transition 
of M = 4, SU (N) supersymmetric Yang-Mills theory on 1 x S 3 was mapped by the 
duality correspondence to the previously discovered Hawking-Page phase transition in 
an anti-de Sitter background [8]. 

With a field theory intuition in mind, a new prediction for the higher dimensional 
black holes (and a phase transition in an infinite volume) was obtained in [4] . Namely, 
using the fact that the chirally symmetric phase of the Klebanov-Strassler (KS) cas- 
cading gauge theories [9, 10] exists only above a certain critical temperature, it was 
proposed in [4] that the dual supergravity backgrounds also would have a regular 
Schwarzschild horizon only above certain horizon temperatures. While there is a good 
understanding of the high temperature thermodynamics of the system (in a chirally 
symmetric phase) [5], the low temperature phase (and the phase transition) requires 
substantial numerical work and is not yet understood. Thus, strictly speaking, the 
prediction for the new type of black holes of [4] has not yet been verified in the KS 
model. Nevertheless, black holes that exhibit a phase transition predicted in [4] were 
shown to exist [6, 7] in a different (but closely related) model — the supergravity dual 
to pure Af — 1, SU(N) SYM theory [11]. Unfortunately, the latter black holes were 
shown to be thermodynamically unstable [12, 7] 1 , and thus it is not clear whether the 
phase transition, while mathematically allowed, is actually physically occurring. 

In this paper we discuss yet another system which we argue undergoes a finite 

1 These black holes have a negative specific heat. 
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temperature phase transition, namely the M = 2* theory, or in other words M = 4, 
SU(N) SYM softly broken by a hypermultiplet mass term to M = 2 gauge theory. 
The supergravity dual to this gauge theory was constructed by Pilch and Warner (PW) 
in [13], and the precise duality map between the gauge theory and the supergravity was 
explained in [14,15]. Here we consider the M = 2* gauge theory at finite temperature, 
realized as a non extremal deformation of the PW flow. The argument for the existence 
of a phase transition in this system is quite simple. For the high temperature phase, 
the mass deformation is irrelevant and we expect the standard near-extremal D3 brane 
thermodynamics. In this phase (which we refer to as the "black hole phase", BH) all 
the thermodynamics quantities scale as iV 2 ; for example, the entropy goes as oc 
N 2 VT 3 . On the other hand, the low temperature phase is rather different. Imagine first 
completely turning off the temperature. In this case the M = 2* gauge theory is exactly 
soluble [16]. Its low energy effective description (valid well below the hypermultiplet 
mass scale) is given in terms of M = 2, {/(l)-^ -1 gauge theory. This gauge theory has 
on the order of iV free degrees of freedom. We expect that this low energy effective 
description is still valid for temperatures much lower than the hypermultiplet mass. 
Thus in the low temperature phase (which we denote the "finite temperature Coulomb 
phase" or the Pilch- Warner phase, PW) the thermodynamic quantities of the M = 2* 
gauge theory would scale as the first power of N, so that, e.g., for the entropy, Spw oc 
NVT 3 . This iV 2 versus iV scaling of degrees of freedom suggests the presence of a 
phase transition 2 between the BH and PW phases, each of which in principle exists 
at all temperatures. 

The physical picture we have in mind for the phase transition is as follows. Start 

with a high temperature BH phase. In this phase the free energy will start negative 

at very high temperatures (as for the black D3 branes), and would gradually increase 

as the temperature is lowered. We expect that at some T = T c the free energy would 

become zero, Fbh{T c ) = 0, and for T < T c it would become positive. In the large N 

limit, the free energy of the PW phase is down by 0(1/ N), and may be taken to be zero. 

Thus the PW phase should be thermodynamically favorable for low temperatures, 

T < T c . We expect that both the PW and the BH phases are thermodynamically 

stable. Moreover, there should be a well defined high temperature expansion of the 

2 We expect a phase transition in the strict N — > oo limit. It is likely that at finite N the phase 
transition is replaced by a crossover regime. We thank Arkady Vainshtein for a useful discussion on 
this point. 
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BH phase 3 . 

In the next section we summarize our results. In section 3 we recall the salient 
features of the Af = 2* gauge theory, and discuss the expectation for its finite tempera- 
ture deformation. We then review the PW Af = 2* renormalization group (RG) flow in 
five dimensions and discuss its non-extremal deformation. The exact ten dimensional 
lift of the "temperature deformed" five dimensional PW flow is constructed in section 

4. In section 5 we analytically determine the leading in m/T correction to the near- 
extremal D3 brane geometry induced by the Af = 2 hypermultiplet mass m. Finally, 
in section 6 we discuss the thermodynamics and the signature of the phase transition. 
First, we explain the computation of the free energy and the energy (mass) of the 
deformed supergravity backgrounds, and compute the difference of the free energies of 
the BH and the PW phases. This computation is valid for arbitrary values of m/T, 
and thus can be used to study the phase transition. Then, using the results of section 

5, we compute the leading correction (in the high temperature phase) to the black D3 
brane thermodynamics. We find that the first law of thermodynamics applied to the 
high temperature phase is violated. We speculate on the relevance of the (induced) 
chemical potential for the resolution of this paradox. We end with some comments on 
the numerical verification of the phase transition. 

Before proceeding, we would like to comment on the study of the thermodynamics 

of the closely related Af = 1*, SU(N) gauge theory [18]. As its name suggests, the 

Af = 1* gauge theory is Af = 4, SU(N) SYM softly broken by a chiral multiplet 

mass term to Af = 1. Unlike the Af = 2* theory, however, the dual supergravity 

background to this model [19] (PS) is known only in the probe approximation. The 

study of thermodynamics in [18] was also done in the probe approximation. Only 

the entropy was computed in the high temperature regime of the non-extremal PS 

background; while the free energy and the energy were not computed independently, 

they were obtained by enforcing the first law of thermodynamics. As we will see from 

the thermodynamics of the Af = 2* model discussed here, the computation of the 

entropy alone does not allow us to reproduce the free energy — it appears one needs 

to compute the induced chemical potential as well. It would be of interest to repeat 

the Af = 2* analysis presented here to the Af = 1* model. But first, the exact extremal 

geometry of the Af = 1* theory has to be understood. 

3 The high temperature expansion of the KS model developed in [5] is ill defined in the ultraviolet. 
This is related to the unusual UV properties of the cascading gauge theories; for a review see [17]. 
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2 Summary of results and outlook 



As the bulk of the paper is rather technical, we highlight our main results in this sec- 
tion. 

■ We first observe that the five dimensional gauged supergravity flow of J\f = 2* 
PW [13] can be deformed to yield a non-extremal black hole geometry with regular 
horizon. The consistency of the D — 5, H — 8 gauged supergravity truncation then 
implies that this black hole solution can be uplifted to the full ten dimensional solution 
of type IIB supergravity. We explicitly verify that this is indeed so. This non-extremal 
deformation is interpreted as the supergravity dual to the finite temperature J\f = 2* 
SU(N) gauge theory in the deconfined phase, which we refer to as the BH phase. 

■ We show that there is a three parameter family of five dimensional black holes admit- 
ting regular horizons. These three parameters are the temperature T and the (generi- 
cally different) masses of the bosonic m& and fermionic m/ components of the M = 2 
hypermultiplet. All regular horizon non-extremal solutions asymptote to AdS§, which 
is consistent with the gauge theory expectation that both the temperature and the 
mass deformations should be irrelevant in the ultraviolet of the gauge theory. Asymp- 
totic M = 2 supersymmetry of the extremal PW geometry imposes a constraint on the 
leading nontrivial asymptotics of the two five dimensional supergravity scalars in the 
non-extremal deformation. The latter reduces the number of independent parameters 
of the regular horizon solution to two: one related to the temperature, and the other 
to the M = 2 hypermultiplet mass m = m\> — nif. 

■ In the high temperature limit, ol\ = (m^/T) 2 <C 1, a 2 = m f/T <C 1, the five di- 
mensional black hole solution is a small deformation of the finite temperature AdS$ 
geometry, representing the S 5 reduction of the throat region of the near extremal D3 
branes. We analytically determine the leading correction in ai of the near extremal 
AdS$ geometry. As expected from gauge theory arguments, asymptotic M = 2 super- 
symmetry sets «i ~ (a 2 ) 2 . 

■ After constructing the black hole solution, we turn to the study of thermodynamics. 
We discuss the computation of the free energy F, the entropy S, and the energy E of the 
non-extremal deformation of the PW flow. The entropy is just the Bekenstein-Hawking 
entropy of the horizon, and is determined from the infrared data of the geometry. The 
free energy, or more precisely FT, is the Euclidean gravitational action, and the energy 
E is the conserved ADM mass of the geometry. Note that computing both E and F 
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requires the knowledge of the IR and the UV data of the solution. Furthermore, we 
verify that F = E — TS is identically satisfied in the supergravity. 

■ While the computation of the entropy is straightforward, both the free energy and 
the energy diverges, and requires regularization. Following [20], we compute F and 
E with respect to a reference geometry which we take to be the supersymmetric PW 
flow with periodically identified (Euclidean) time direction with periodicity equal to 
the inverse horizon temperature. We call this geometry the U PW phase". The pre- 
scription of [20] requires the introduction of a boundary cutoff and the matching of 
induced geometries (and matter fields) for the background at hand and the reference 
one "up to sufficiently high order" [20]. We apply the "minimal subtraction" prescrip- 
tion for matching, where only the leading asymptotics of the induced geometries and 
the matter fields are matched. This prescription gives the correct answers for simple 
black hole geometries such as the Schwarzschild-anti-de Sitter solution. It also works 
in more complicated cases such as the nonabelian black hole solutions of [7]. 

■ Using this minimal subtraction prescription for the free energy and the energy, we 
find an explicit analytical expression for F BH — F PW (or E BH — E PW ) is terms of the 
coefficients of the subleading ultraviolet asymptotics of the five dimensional scalars in- 
ducing the Pilch- Warner flow [13]. An added bonus of using the PW background (with 
appropriately compactified Euclidean time direction) as the reference one in F and E 
regularization is the fact that the purported phase transition between the high temper- 
ature BH phase and the low temperature PW phase arises when Ap^, = F BH — F PW 
changes sign. 

■ Using the high temperature expansion (corresponding to deformations from the 
non-extremal AdS^ geometry), and working to leading order in m/T, we analyti- 
cally compute the corresponding deformations of the thermodynamic quantities. While 
F = E — TS continues to be satisfied after regulation, we however find that TdS ^ dE. 
We have verified our prediction for the leading correction to the free energy numeri- 
cally. This indirectly confirms the violation of the first law of thermodynamics. 

■ Though we have been unable to find a satisfactory explanation for this apparent con- 
tradiction with the first law of thermodynamics for the high temperature phase of the 
J\f = 2* flow, we point out that this paradox could be resolved once we include a certain 
chemical potential induced by the fermionic mass term of the M = 2 hypermultiplet. 

Perhaps the most intriguing conclusion we have reached is the fact that the proper 
interpretation of a finite temperature deformation of the Pilch- Warner geometry ap- 
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pears to require the introduction of a nonvanishing chemical potential dual to the 
sources that are turned on in the UV. This induced chemical potential follows from the 
conjecture that the string theory partition function in the gauge/string theory corre- 
spondence is dual to the gauge theory grand canonical partition function. In section 
6.3 we outline the general arguments leading to such a statement. It would be inter- 
esting to verify this in a more general setting, e.g., by studying the finite temperature 
deformations of generic holographic renormalization group flows as in [21,22]. 

The original motivation for the study of the Af = 2* thermodynamics presented 
here was to study and confirm the phase transition between the BH and the PW 
phases. However, our analysis for this question is as yet inconclusive. While the high 
temperature expansion is suggestive that such a phase transition occurs, additional 
analytical or numerical work is required to extrapolate to the region m ~ T where 
solid evidence of the transition would be obtained. We hope to report on these results 
in a separate publication. Finally, it is interesting to understand the "hydrodynamic 
description" of the BH phase of the Af = 2* gauge theory along the lines of [23]. 

3 J\f = 2* RG flow and its non-extremal deformation in five 
dimensions 

3.1 The gauge theory picture 

In the language of four-dimensional Af = 1 supersymmetry, the mass deformed Af = 4 
SU (N) Yang-Mills theory (Af = 2*) in R 3,1 consists of a vector multiplet V, an adjoint 
chiral superfield $ related by Af = 2 supersymmetry to the gauge field, and two 
additional adjoint chiral multiplets Q and Q which form an Af = 2 hypermultiplet. 
In addition to the usual gauge- invariant kinetic terms for these fields 4 , the theory has 
additional interactions and a hypermultiplet mass term given by the superpotential 

W = ~T~ Tr([Q, Q]$) + -4r-(Tr<2 2 + TrQ 2 ) . (3.1) 

9ym 9ym 

When m = the gauge theory is superconformal with g YM characterizing an exactly 
marginal deformation. The theory has a classical 3(iV — 1) complex dimensional mod- 
uli space, which is protected by supersymmetry against (non)-perturbative quantum 
corrections. 

4 The classical Kahler potential is normalized according to (2/(?y M ) Tr[$$ + QQ + QQ]. 
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When m ^ 0, the M = 4 supersymmetry is softly broken to M = 2. This mass 
deformation lifts the {Q, Q} hypermultiplet moduli directions, leaving the (N — 1) 
complex dimensional Coulomb branch of the M = 2, SU(N) Yang-Mills theory, pa- 
rameterized by expectation values of the adjoint scalar 

$ = diag(ai,a 2 , • • • ,a N ) , ^a; = 0, (3.2) 

i 

in the Cartan subalgebra of the gauge group. For generic values of the moduli a^, 
the gauge symmetry is broken to that of the Cartan subalgebra [/(l)^ -1 , up to the 
permutation of individual U(l) factors. Additionally, the superpotential (3.1) induces 
the RG flow of the gauge coupling. While from the gauge theory perspective it is 
straightforward to study this M = 2* theory at any point on the Coulomb branch [16], 
the PW supergravity flow [13] corresponds to a particular Coulomb branch vacuum. 
More specifically, matching the probe computation in gauge theory and the dual PW 
supergravity flow, it was argued in [14] that the appropriate Coulomb branch vacuum 
corresponds to a linear distribution of the vevs (3.2) as 

aie[-a ,a ], a 2 = — — , (3.3) 

7T 



2 / f a ° 

P( a ) — — — \/ a o~ fl2 > / dap(a)=N. (3.4) 
m 9ym V J-an 



with (continuous in the large N limit) linear number density 

"do 

m "9YM V J-a 

Unfortunately, the extension of the N = 2* gauge/gravity correspondence of [13-15] 
for vacua other than (3.4) is not known. 

In [14, 15] the dynamics of the gauge theory on the D3 brane probe in the PW 
background was studied in detail. It was shown in [14] that the probe has a one complex 
dimensional moduli space, with bulk induced metric precisely equal to the metric on the 
appropriate one complex dimensional sub manifold of the M = 2*, SU(N + 1) Donagi- 
Witten theory Coulomb branch. This one dimensional submanifold is parameterized 
by the expectation value u of the U(l) complex scalar on the Coulomb branch of 
the theory where SU(N + 1) — > U(l) x SU(N)pw- Here the pw subscript denotes 
that the SU(N) factor is in the Pilch- Warner vacuum (3.4). Whenever u coincides 
with any of the otj of the PW vacuum, the moduli space metric diverges, signaling the 
appearance of additional massless states. An identical divergence is observed [14, 15] 
for the probe D3-brane at the enhangon singularity of the PW background. Away 
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from the singularity locus, u = a e [— a ,a Q ], the gauge theory computation of the 
probe moduli space metric is 1-loop exact. This is due to the suppression of instanton 
corrections in the large iV limit [14, 24] of M = 2 gauge theories. 

Consider now M = 2* gauge theory at finite temperature T. Turning on a mass 
m for the hypermultiplet sets a strong coupling scale A oc m. We expect to find two 
different phases of this gauge theory, depending on whether T > A or T < A. In the 
former case the effect of the mass deformation is negligible, and we expect to recover 
the M = 4 thermodynamics. In particular, conformal invariance dictates that the free 
energy scales like T 4 , with a prefactor of iV 2 indicative of the scaling of the number of 
degrees of freedom. At weak 't Hooft coupling the familiar result reads [25] 

Fsym = -yiVVT 4 . (3.5) 

By symmetry arguments, we expect the corrections to the free energy (3.5) due to 
the mass deformation (3.1) to be of order O (m 2 /T 2 ). The A^ 2 scaling of the thermo- 
dynamic quantities naturally occurs in the ten-dimensional black holes describing the 
non-extremal deformation of the dual supergravity backgrounds. For this reason we 
will call the high temperature phase of M = 2* SYM the BH or black hole phase. 

In the other limit, we expect qualitatively different physics in the low temperature 
phase of M = 2* SYM. Ignoring T in the first approximation, the low energy effective 
description of the Af = 2* theory is given by free [/(l)^ 1 SYM as explained above. 
This effective description breaks down at scales of order the strong coupling scale, i.e. 
m, but is appropriate as we turn on the temperature provided T <^ m. The number 
of (free) degrees of freedom of this effective low energy description scales like N, which 
must be reflected in the scaling of the thermodynamics quantities such as the entropy, 
S oc NVT 3 . We denote this phase the "finite temperature Coulomb phase", or the 
PW phase. Notice that 

FpW ~ - (3.6) 



Fbh AT' 

and thus vanishes in the large- A" limit. Qualitatively, we expect the free energy of the 
J\f = 2* SYM to behave as in Fig. 1. 

Given two possible phases, BH and PW, the most favorable one is that with lowest 
free energy. Thus the signature of a phase transition in the M = 2* theory would be 
a change in sign in the difference Ap^ = Fbh ~ Fpw- This suggests that, in order to 
have a phase transition, the sign of the 0(m 2 /T 2 ) correction to (3.5) must be positive 
so that, when the temperature is lowered from an initial high temperature phase, Fbh 
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Figure 1: The expected behavior of free energies of the BH (red curves) and the PW 
(blue curve) phases (Fbh arid Fpw) as a function of T/m in the large N limit. The 
dotted line for the BH phase would imply a second order phase transition, while the 
solid line would imply a first order phase transition. 

will be driven above Fpw at a finite temperature T c . While we have not performed 
this computation perturbatively in the 't Hooft coupling, we have instead analytically 
determined the 0(m 2 /T 2 ) correction for strong 't Hooft coupling in the dual supergrav- 
ity computation. We find that this correction indeed has a positive coefficient (6.48) 
at strong 't Hooft coupling. This provides evidence for a phase transition, although 
further investigation is necessary beyond leading order to substantiate this claim. 

3.2 The PW renormalization group flow 

The gauge theory RG flow induced by the superpotential (3.1) corresponds to a five 
dimensional gauged supergravity flow induced by a pair of scalars, a = In p and x- 
The effective five- dimensional action is 



S 



AttG, 



1 




(3.7) 
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where the potential V is 5 



1 

16 



dW 
da 



with the superpotential 



+ 



1 



dW 

dx 



W = — - - -p 4 cosh(2x) 



(3.8) 



(3.9) 



pr 2' 

Note that we have chosen identical normalizations for a and x to highlight the general 
D — 5, M — 2 features of the flow. The PW geometry [13] has the flow metric 



ds\ = e 2A (-dt 2 + dx 2 ) + dr 2 



(3.10) 



Solving the Killing spinor equations for a supersymmetric flow then yields the first 
order equations 



— --W da_\dW_ dx_ldW 
dr 3 dr 4 da ' dr 4 dx 



(3.H) 



It is straightforward to verify that solutions to these flow equations will automatically 
satisfy the scalar and Einstein equations of motion. 



3.2.1 Asymptotics of the PW flow 

This system was solved in [13] by rewriting the equations in terms of x as an indepen- 
dent variable. Given this explicit solution of the flow equations, it is easy to extract 
the UV/IR asymptotics. In the ultraviolet, r — * +oo, we find 

UV: X ^0 + , A^^r. (3.12) 

This corresponds to the scalars approaching the maximally symmetric AdS 5 x S 5 UV 
fixed point of the potential V . In the infrared, r — * 0, we find instead 

IR: p^0 + , x^+oo, A ^-\*- (3-13) 

As will be apparent later, this flow to the IR will be cut off at finite temperature. 

5 We set the 5d gauged supergravity coupling to one. This corresponds to setting the S 5 radius 
L = 2. 
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3.3 The non-extremal PW flow 



We now consider deforming the PW flow by turning on non-extremality in the metric. 
Since the deformed flow breaks supersymmetry, we can no longer appeal to first order 
equations, but must consider the second order equations of motion. The action (3.7) 
yields the Einstein equation 

^R,u, = d^ad v a + d^ X d vX + ^9^V , (3.14) 

and the scalar equations 

ldV l&P , 01r . 

Da= 2^' DX =2dx- (3 - l5) 

For a finite temperature deformation of the flow metric (3.10), we take 

dsj = e 2A (-e 2B dt 2 + dx 2 ) + dr 2 , (3.16) 

where e 2B represents a blackening function. Note that we choose to retain g rr = 1 since 
any non-trivial factor can be absorbed into a redefinition of r. 

Substituting this metric ansatz into the equations of motion, (3.14) and (3.15), we 
find 

= B" + {4A' + B') B' , (3-17) 
\A" + \ B " + (Af + \ (Bf + \a>B> = -hp , 
- A" - \ B " - (Af - \ (B') 2 - l -A>B> = (a') 2 + ( X f + hp . 

Note that, by defining A = A + \B and B = \^3B/A and taking appropriate linear 
combinations, the above equations may be written in the equivalent form 

1 dV 

= a" + 4A'a'-^, 

2 da 

1 dV 

2 OB ' 

3(A') 2 = (a') 2 + (x') 2 + (B') 2 -V, 
3A" = -4 ((a') 2 + { X 'f + {B'Y) , 
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where we have formally introduced dV/dB = 0. When written in this form, it is easy 
to see that the last equation is redundant, and may be obtained by differentiating the 
penultimate equation and substituting back in the scalar equations of motion. Since 
these equations are consistent, we can use the same scalars as in the PW case, even 
when considering deformed flows. 

At this point, a few comments are in order. Firstly, the last equation of (3.18) 
provides a non-extremal generalization of the holographic c-theorem, namely A" < 
where e 4A = \f—g. Secondly, scalars in AdS 5 may be labeled by the representations 
D(E ,0, 0) where E is the lowest energy state, and may be related to the conformal 
dimension, A, of the dual field theory operators. Expansion of V about the UV fixed 
point indicates that a (or p) has E Q = 2, while x has E = 3. The blackening factor 
B may be thought of as a scalar mode with E = 4. Finally, we see that the equation 
for B in (3.17) can be integrated once to obtain 



or equivalently In B' + A A = const. This relation will prove useful below. 
3.3.1 Asymptotics of the finite temperature deformation 

The thermal solutions we are interested in have regular black hole horizons. As a result, 
we may examine the solution to the system of equations (3.17) near the horizon. In 
this case, the behavior of (3.13) is cut off, and the scalars run to fixed values, po an d 
Xoj 011 the horizon. Including the horizon value of A, we see that the nonsingular in 
the IR flows are given by a three parameter family {a, po > 0, Xo}> specifying the near 
horizon (r — > 0) Taylor series expansions 



\nB' + 4A + B = const, 



(3.19) 




(3.20) 
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Here, 5 = 5(p ,Xo) should be adjusted so that e B — > 1_ as r — > +00. The first 
non-trivial terms in the series expansions (3.20) are 



(3.21) 



cT 2 (I1 — J2 Po 4 + I Po cosh(2xo) - ^ Po sinh 2 (2xo) , 
cT 2 61 = -§ po 4 - § po cosh(2xo) + ^Po sinh 2 (2xo) , 
Pi = i Po 3 - k Po cosh(2xo) + ^ Po sinh 2 (2xo) , 
5- 2 xi = -| Po sinh(2xo) + ^ Po sinh(4x ) • 

4 The ten-dimensional solutions 

In this section, we lift the deformed PW flow to ten dimensions. However before doing 
so we establish our conventions and review some of the pertinent aspects of the lifting 
procedure. 

4.1 Type IIB supergravity equations of motion 

We use a mostly positive convention for the signature ( — h • • • +) and take ei...io — + 1. 
The type bosonic IIB equations consist of the following [26]: 
• The Einstein equations: 

Rmn = Tm N + Tm N + Tm N , (4.1) 

where the energy momentum tensors of the dilaton/axion field, B, the three index 
antisymmetric tensor field, F( 3 ), and the self-dual five-index tensor field, F( 5 ), are given 
by 

(4.2) 

Tmn = g(G PQ mG* pqn + G* PQ m Gpqn - q9mnG pqr G* pqr ) , (4.3) 

and 

= \f PQRS mF pqrsn . (4.4) 
In the unitary gauge, B is a complex scalar field, and 

P M = f 2 d M B , Q M = f Im (Bd M B*) , (4.5) 

where 

f = (1 - BB*) 1 / 2 ' (4 ' 6) 
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while the antisymmetric tensor field G( 3 ) is given by 

G ( 3) = /(F ( 3)-£F ( * 3) ). (4.7) 

• The Maxwell equations: 

(V p — iQ P )G M NP = P P G* MNP — - i F MNP q R G pqr . (4.8) 

• The dilaton equation: 

(V M - 2tQ M )P M = -±G PQR G PQR . (4.9) 

• The self-dual equation: 

F (5) =*F {5) . (4.10) 

In addition, F( 3 ) and F( 5 ) satisfy Bianchi identities which follow from the definition of 
the field strengths in terms of their potentials: 

F( 3 ) = dA {2 ) , 

1 (4.H) 
F {5) = dA (A) - -Im(A (2) A F ( * 3) ) . 

For the ten-dimensional uplift of the RG flows in the five-dimensional gauged super- 
gravity, the metric ansatz and the dilaton is basically determined by group theoretical 
properties of the D — 5, M — 8 scalars. Thus they must be the same for both the 
deformed and original PW flows. Specifically, we assume [13] that the D — 10 Einstein 
frame metric is 

ds?0 = n*d4 + 4 (J^iMH L-w + p » cos2 e (4- + £±£) + S i„- e d -f) , 

p \ \cX 2 Xi J X 2 J 

(4.12) 

where ds 2 is either the original PW flow metric (3.10) or its deformations (3.16), and 
c = cosh(2x). The warp factor is given by 



, 2 (cX.X,)^ 



n 2 = , (4.i3) 

p 



and the two functions X { are defined by 

Xi(r, 9) = cos 2 9 + p{rf cosh(2x(r)) sin 2 9 



(4.14) 

X 2 (r, 9) = cosh(2x(r)) cos 2 9 + p(r) 6 sin 2 9 . 



15 



As usual, <7j are the SU(2) left-invariant forms normalized so that doi = 2ajAa^. Note 
that we perform all computations in the natural orthonormal frame given by 

e 1 oc dt, e 2 oc dr, e 3 oc dxi, e 4 oc dx 2 , e 5 oc dx%, 

(4.15) 

e 6 oc d9, e 7 oc <Ti, e 8 oc a 2) e 9 oc <t 3 , e 10 oc g?0 , 
Turning now to the matter fields, for the dilaton/axion we have 

'-s((f)"*($r> '^((tr-(trv 

(4.16) 

The consistent truncation ansatz does not specify the 3-form nor 5-form fluxes. As 
in [13], for the 2-form potential we assume the most general ansatz allowed by the 
global symmetries of the background 

A {2) = e i0 (ai(r,0) d9 Aa 1 + a 2 (r,9) a 2 Aa 3 + a 3 (r,9) a x Arf0 + a 4 (r, 9) d9 Ad(f) , (4.17) 

where aj(r, 9) are arbitrary complex functions. For the 5-form flux we assume 

F 5 = T + *T, J- — dt A vol M a A du , (4. 18) 

where u(r, 9) is an arbitrary function. As in the PW case, examination of the Einstein 
equations reveals that 2-form potential functions a, have the following properties: = 
0; ai, a 2 are pure imaginary, and 03 is real. 

4.2 Lift of the near extremal deformation 

The verification of the uplifted solution proceeds exactly as for the M = 2* flow defor- 
mations discussed in [27]. Thus we present only the results. We find 

ai = — i 4 tanh(2%) cos#, 

. p 6 sinh(2x) . Q 2q 
a 2 = iA — sin cos 9, (4.19) 

sinh(2x) . o 9, 



a 3 = —4 — sin 9 cos 9 , 

Ao 



and 



^ = l - e iA+B 1 (-p 12 sinh 2 (2 X ) sin 2 9 + 2p 6 cosh(2 X ) (1 + sin 2 9) + 2 cos 2 ^ . 

(4.20) 



16 



We have explicitly verified that by supplementing the metric and the dilaton/axion 
ansatz of the previous section with (4.19), (4.20) and the five- dimensional flow equa- 
tions (3.17), all the equations of ten-dimensional type IIB supergravity are satisfied. 

5 High temperature expansion 

Having examined the system of equations governing the non-extremal flow, (3.17), we 
now turn to the construction of solutions. At finite temperatures, we find it convenient 
to parametrize the flow not in terms of the radial coordinate r, but rather in terms of 
the blackening function e B . To do so, we introduce a new coordinate 

y = e B , ye [0,1], (5.1) 

with y = being the horizon and y — > 1_ the UV asymptotic limit. The standard near- 
extremal D3 brane solution is realized when the bosonic and fermionic masses of the 
J\f = 2 hypermultiplet components are turned off, corresponding to the supergravity 
scalars p and x sitting at the UV fixed point. The near-extremal D3 brane solution 
has the form 

A(y) = & - Iln(l - y 2 ) , 

P(y) = i, ^ 
x(y) = o, 

where a is an integration constant which determines the BH temperature according to 

T = —e & . (5.3) 
2vr V ; 

We recall that at zero temperature the PW flow involves the scalars p and x running 
away from the UV fixed point as one flows to the IR. At a regular horizon, the scalars 
attain fixed values po and xo- Hence we now seek a solution to (3.17) satisfying the 
conditions 

(P,X)^(!,0) as y->l_, (p,x) -»• (f>o,Xo) as 3/ . (5.4) 

Several flows satisfying these boundary conditions are displayed in Fig. 2. It is clear 
from the figure that the flows proceed further into the IR as the temperature is lowered. 
While we have been unable to find an exact analytical solution, it is nevertheless 
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a 

Figure 2: Finite temperature PW flows away from the UV stable fixed point (a, x) = 
(0, 0). The endpoints of the flows correspond to the fixed values of the scalars at the 
horizon. The dashed line gives the supersymmetric zero-temperature flow. 




possible to develop a consistent (uniformly convergent) perturbative approximation in 
the high temperature phase. On the gauge theory side, this corresponds to a power 
series expansion in oti oc (m^) 2 <C 1 and a 2 oc rrif <C 1, where nif, and rrif are masses 
of the bosonic and fermionic components of the M = 2 hypermultiplet measured with 
respect to the string scale. In what follows, we solve for the leading order deformation 
in «i and a-i. Specifically, we seek a solution of (3.17) in the form 

A(y) =a-~ ln(l - y 2 ) + afA^y) + a 2 2 A 2 (y) , 

p{y) = l + a lPl {y), ( 5 - 5 ) 

X(y) = "2X2(2/) • 

Substituting this ansatz into (3.17), and working to first non-trivial order in Oj, we 
find the linearized scalar equations 

= (l-y 2 ) 2 (yp'J + yp!, 

= (1 -2/ 2 ) 2 {yx'2)' + \vx2, (5 ' 6) 

as well as the equations governing the back-reaction on the metric 
= y(l - y 2 ) A>[ - (1 + 3y 2 ) A\ + %(1 - y 2 ) (p[) 2 , 

= y(l- y 2 ) K - (1 + 3y 2 ) A' 2 + iy(i - y 2 ) ( X ' 2 ) 2 . (5 ' T) 



The scalar equations, (5.6), may equally well be obtained from the linearization 
of (3.15). Note that for an arbitrary scalar $(?/) of mass m, its equation of motion, 
□ $ = m 2 Q, in the background (5.2) has the form 

$" + 1$' = ( mL ) 2 $ . (5.8) 
V 4(1-1/2)2 ^ 

This may be readily solved in terms of hypergeometric functions. Although there are 
generally two linearly independent solutions, only one combination is regular at the 
horizon, y — » 0. Defining E = 2 ± a/4 + (mL) 2 , the regular solution has the form 

$ = (l-y 2 ) E « / SF i e i E ,lE ,l;y*). (5.9) 

As a result, for the p and x scalars, we obtain 

pi = (i -v 2 ) 112 2 *\(y,i;y 2 ), 
x 2 = (i-y 2 ) 3/4 2 ^i(I,| ? i^ 2 )> 

where, without loss of generality, we assumed the horizon boundary conditions 

Pi 



= 1, X2 

y=0 



= 1. (5.11) 

y=0 

Note that as y — > 1_, the perturbations pi, %2 vanish. This is readily seen by rewriting 
the solution, (5.9), as 

$ = rn"iF°| ( 1 - ^Vi^o, |^o, l - y 2 ) 

r(1 "^ (5.12) 

+ W^ i 1 " ^ 2 ) (4 - £o)/4 2 ^i(l - i^o, 1 - l^o, 2 - §£ ; 1 - y 2 ) , 

which is valid provided 6 E 7^ 2. Expanding for y — > 1_ yields the expected boundary 
behavior, $ ~ (1 — y) A+ ^ and (1 — y) A ~/ 4 , where the conformal dimensions A± corre- 
spond to Eq and 4 — £0- 111 this case, however, instead of being independent, the A + 
and A_ modes are related by the condition of horizon regularity. 

Turning to the leading order gravitational back-reaction, we see that (5.7) can be 
solved by quadratures 



(1 — z 2 ) 2 \ In \ dx J X 
_ v (5.13) 

4={ ^fA^f^^'' 1 - 111 " 



3 Jo (1 — z 2 ) 2 \ J \ dx J x 
6 For pi, which has Eq = 2, the behavior as y — > 1_ picks up a log, namely pi ~ (1 — y) 1 ^ 2 and 
(1 - y) 1 ^ ln(1 _ y) . 
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where 7,, £j are four integration constants. For a generic choice of 7$, we find A i \ y _> 1 _ oc 
(1 — thus to recover the proper asymptotics in the UV AdS§ x S 5 geometry, these 
constants must be fine tuned: 

8-7T 2 8 - 37T . 

The other two integration constants, £j, can be absorbed in a redefinition of a. In 
fact, we show below that the physical quantities are independent. Note that at the 
horizon we have the behavior 

= 5 2 ~ 1 + 4aj (1 + 2 7 i) + + • (5.15) 



dy_ 

This allows us to compute the BH temperature 



3 



2vr Ur 



27 6 ° 5 ~ 2^ f 1 + (2 + 6 + 47l) + ^ (3 + 66 + §72) 



(5.16) 

While the first relation is valid for arbitrary temperature and masses, in the second 
one we have kept only the leading terms in CKj. 

Using the explicit lifting of the metric, (4.12), we may compute the area of the BH 
horizon 



A _ T/ P 3A 

•^horizon v 3 c 



= V 3 e 3a 

v-+o+ (5.17) 
~ V 3 e 3& 2 5 vol 55 (1 + 3^« 2 + 3&c$ , 

where V3 is the 3-dimensional volume and V0I55 is the volume of the unit S 5 . Again, the 
first relation in (5.17) is exact for all temperatures/masses. The Bekenstein-Hawking 
entropy density is 7 

S BB = ^ = \ M Q- e «) 3 ~\ M Q- e*) 3 (1 + % «? + 3^) 

(5.18) 



6 Thermodynamics and the signature of the phase transition 

In this section we discuss the thermodynamic properties of the N = 2* theory. As we 

have explained above, physically we expect a phase transition between the deconfining 

7 We have used the standard relations 16ttGn — (2n) 7 gglf, 4:Trg s Nl 4 — L 4 , and the fact that we set 
L = 2. 
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phase (at high temperature) and the finite temperature Coulomb phase (at low tem- 
perature). The high temperature phase is realized by the BH geometry, represented by 
the solution to (3.17) with boundary conditions (3.20). The low temperature phase is 
the 7 = (Euclidean) PW geometry [13] with periodically identified (Euclidean) time 
direction t E ~ t E + 1/T. 

We begin by considering the standard definition (and regularization) of the free 
energy and the energy of the finite temperature deformed PW background. Specifically, 
we identify the Helmholtz free energy F with the combination T I E enom , where J^ nom 
is the renormalized Euclidean gravitational action, and the dual gauge theory energy 
E with the ADM mass of the finite temperature deformed PW geometry. Also, we 
identify the gauge theory entropy with the Bekenstein-Hawking entropy of the deformed 
PW background. We show that, with such identifications, we identically satisfy the 
thermodynamic relation F = E—TS. This extraction of the thermodynamic quantities 
are valid for arbitrary values of mass and temperature. 

To proceed, we note that supersymmetry of the PW background relates the a\ and 
«2 coefficients of the leading nontrivial asymptotic behavior of the five-dimensional 
scalars p and x- This allows us to parametrize the thermal phase of PW by the 
single quantity m/T. In the high temperature phase, we analytically compute the 
leading correction to the near extremal D3 brane thermodynamics due to the PW 
mass deformation. The sign of this correction to the black D3 branes is consistent with 
our claim for a phase transition. 

We find that, in the high temperature phase, our extracted free energy no longer 
satisfies dF = —SdT, thus apparently violating the first law of thermodynamics. We 
provide a possible resolution to this puzzle in terms of a generalized chemical potential 
induced by the PW mass deformation. However a full understanding of the thermody- 
namics requires additional investigation. 

6.1 The regularized free energy and the energy 

We recall that the free energy F, the energy E, and the entropy S of the a system is 
related by the well known expression 

F = E-TS. (6.1) 

Here, F, E and S are well defined quantities in a weakly coupled gauge theory, and 
physically should remain well defined (finite) at large 't Hooft coupling. It is known, 
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however, that in the dual supergravity both the free energy and energy densities are 
divergent, and thus need to be properly renormalized before they can yield a finite 
answer. A standard regulating procedure is to compute such quantities by compar- 
ison with a 'reference' supergravity background having the same asymptotics. This 
comparison is often ad hoc; in particular, the matching of the two geometries at hand 
at the regularization boundary remains slightly ambiguous. In our case, however, we 
have a physically well motivated reference geometry, namely the PW background with 
periodic Euclidean time direction of appropriate size. 

Using the PW background as reference, strictly speaking we will not be computing 
F BH , E BH and S BB directly [where the subscript BB relates to the non extremal 
deformation (3.16)-(3.20)], but rather the differences 



AF = F BH — F PW , AE = E BH — E PW , AS = S B h — S PW . (6.2) 



In practice, we expect these quantities to be dominated by their BH values, AF = 
F BH , AE = E BH , AS = S BH . This is clearly the case at weak 't Hooft coupling 
since the thermodynamic quantities in the finite temperature Coulomb phase are 1/N 
down compared to the corresponding quantities in the deconfined phase, and thus the 
former are essentially zero in the large N limit. Experience with other examples of the 
gauge/gravity correspondence suggests that going to strong 't Hooft coupling would 
typically modify the prefactor, but not the large iV scaling of the free energy, the 
energy and the entropy. Note that the choice of PW background as a reference one 
is particularly convenient when exploring the phase transition, as a phase transition 
implies going through a zero in AF in (6.2) as one changes the temperature. 

Before proceeding to the BH solution, we recall the asymptotics of the reference 
PW geometry. In [13] the solution of the supersymmetric 7 = flow equations is given 
in terms of x as the flow coordinate: 



The single integration constant k in (6.3) is related to the hypermultiplet mass m in 



sinh(2x) ' 
p 6 = cosh(2x) + sinh 2 (2x) In 



sinh(x) 
cosh(x) ' 



(6.3) 



(3.1) by [14] 



k = mL = 2m . 



(6.4) 
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As indicated in (3.12) and (3.13), the scalars (p,x) start at their UV fixed point 
(1_, + ), and flow toward (0 + , +00) in the IR. Using the the flow equations, (3.11), we 
find the IR asymptotics 

1 !, 4 
Inp ^ _ x+ -In-, 

3 6 3 

e A ~2A;(4/3) 1/3 e ^ /3 , (6.5) 

^~(4/3) 2 / 3 e 2 */ 3 . 
ar 

For matching, we are more interested in the UV behavior. To develop the asymptotics 
at the boundary, we introduce 

x = e- r ' 2 . (6.6) 

We find in the UV 

X ~ kx [1 + k 2 x 2 (§ + § \n(kx)) + k A x A (-£ + f ln(fcr) + f \n 2 {kx))] , 

p ~ 1 + k 2 x 2 (§ + § ln(ifex)) + A; 4 x 4 + 2 ln(fcr) + § ln 2 (A;x)) , 

(6.7) 

A ~ - ln(2x) - |A; 2 x 2 - A; 4 £ 4 (| + f ln(fcr) + | ln 2 (A;x)) , 



cL4 

cir 2 



- ± + |A; 2 £ 2 + A; 4 x 4 (l + | ln(kx) + | ln 2 (fcx)) . 



Note that, as will be evident later, we need to keep terms up to 0(x 4 ) in the expansion. 

Turning now to the BW geometry, the general solution of (3.17) which is smooth 
in the IR (e B — > 0) has three integration constants, {«, Xo, Po} ; which are related to 
temperature and masses of the M = 2 hypermultiplet components, (3.20), (3.21). The 
most general solution of (3.17) in the UV (x — > 0+) has altogether five parameters, 
{£> PiO) Xo> Xio}- Three of them are related to the temperature and the masses, 
while the other two are uniquely determined from the requirement of having a regular 
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horizon, (3.21). In any case, we have a three parameter BH solution 
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B ~ - (3 x 4 [1 + lx 2 x 2 + x 4 ^^ - \ PllPw + ^Xo + 2p? + XoX 
+ \nx (-§p n + |xo + 4p n pio) + 2/3?! In 2 x)] , 

X ~Xo^ [1 + x 2 (xio + fxo lna; ) + ^ 4 (f Pn ~ fPuPio - §Xo ~ §XoPn + 2 XoPio 
+ 5p 2 + 2xoXio + lnx (-f Pn + lOpnPio + §Xo + 2XoPn) + 5p n ln 2 x)] , 

p ~ 1 + x 2 (p 10 + pn lux) + x 4 (-2p n p 10 + fpn + §p 2 + yXoPio - fXoPn + |Xo 
+ In a; (3pnpi + f XoPn ~ 2 Pii) + §Pn ln2 > 

A ~£ - In x - Ixo^ 2 + ^ 4 (i/5 + |xo - IXoXio - |Pn - P? 
- lux (§Xo + 2pnpio) - Pn In 2 x) , 

(6.8) 

Here we have introduced an additional integration constant j3 which, however, can be 
absorbed at the expense of shifting the position of the horizon in the radial coordinate 
r (or alternatively by rescaling x). For this reason, (5 should not be considered an 
independent parameter of the solution. Also, we find 

dA 

-fo ~ \ + l^W + ^ 4 + 2 Plo + PnPio + \p\i + XoXio + \xt ^ ^ 
+ ln x ( P 2 U + fxo + 4pnPio) + 2p 2 n In 2 x) . 

In (6.8) and (6.9), x = x§e~ r l 2 with x an arbitrary constant. 

The free energy, F, of the gravitational action can be obtained from the (Euclidean) 
action I E according to 

F = T I E = -L e a 5 I E [a, Xo, Po] , (6. 10) 

Ztt 

where T is the temperature. As usual, I E is divergent and should be properly regular- 
ized. As explained above, our approach is to regulate the free energy by subtraction, 
AF = TAI E , where 

A/ E [a,xo,Po] = lim \ r BH [a, Xo, Po] - F PW [T,k] \ . (6.11) 
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The regularized action F E consists of both volume and surface terms 



jr _jr _,_ jr 
1 E — J bulk ^ J surf 



l — fdr 1^ d^ v 7 ^ {-\R E + {daf + {d X f + V) (612) 



1 



where G5 is the five dimensional Newton's constant 

G ^2^ = W^ (6 - 13) 

g E is the Euclidean version of the metric (3.16), is a unit vector orthogonal to the 
four-dimensional boundary dM.$, and h® v is the induced metric on dM.^ 

K v =g* v + n^n u . (6.14) 

In (6.11) we have assumed that the boundary dM.*> is defined at fixed r [in the coordi- 
nates (3.16)], which we will take to infinity at the end of the calculations. In this case, 
the unit normal vector is 11^ = 6^. 

Consider first the bulk contribution in (6.12). Because of local diffeomorphism 
invariance, the on-shell value of the action must reduce to a surface integral. This is 
indeed what we find: 8 



(6.15) 



where v is an arbitrary constant parameterizing the constraint (3.19). In what follows, 
we find it convenient to set 

v = 0. (6.16) 

In this case 

(6.17) 



jr _ 1 3A ( p A+B\> 

Ibulk ~8^G~ 5 e [e > 



horizon 

where horizon refers to either the standard black hole horizon location for the "decon- 
fining phase" (BH) analytically continued to Euclidean signature or the IR (x — > +00) 



i We omit the volume integral over the boundary dM$: J gMr d 4 £ = V3/T. 
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of the Euclidean PW solution with periodically identified time direction (PW). Notice 
that at the black hole horizon 



8ttG h 



(e A+B )' 



AA 



horizon, BH 



dy_ 

Or 



— SbhT , 



(6.18) 



horizon, BH 



where S B h is exactly the black hole entropy (5.18), and T is the corresponding black 
hole temperature (5.16). On the other hand, using the IR asymptotics of the PW 
solution, (6.5), we find instead 



„3A 



87rG 5 



— — SpwT , 



(6.19) 



horizon, PW 



which is simply interpreted in terms of the vanishing entropy of the PW phase. 

It should be noted that the black hole horizon is a regular point of the Euclidean 
geometry. Thus it is unusual to find a horizon surface term contribution to the Eu- 
clidean bulk action (6.15). In fact, this contribution is somewhat artificial, and arises 
because of our particular choice of v , (6.16). Indeed, for generic v , we find 



87rG f 



e 3A ( e A +B y + 



V 



AA+B 



4ttG s 



B' 



horizon 



l + 2v 
8ttG 5 



AA+Bgl 



(6.20) 



horizon 



where we have used the fact that e 3A+B (e A )'\ horizon = for both the BH and the PW 
phases. From (6.20) we see that for v = —1/2, the full contribution to (6.15), 
would come from the asymptotic region. Although, strictly speaking, this is the only 
proper value for v, since the full value of the Euclidean action (6.15) is independent of 
v, we nevertheless find it convenient to retain v = 0, as indicated in (6.16). 
For the surface term in (6.12), we find 



V = 

surf 



8irG s 



if 



AA+B 



)' 



(6.21) 



Adding the bulk (6.17) and the surface (6.21) terms together, we find 



,3A 



8ttG s 



{e A+B )' 



AA+B 



horizon 



A' 



(6.22) 



We have shown above that the first term in (6.22) is simply the combination —ST where 
S is the entropy density. If the standard relation (6.1) is realized in the supergravity 
(and it must be so), then the other term in (6.22) must be the regularized energy 
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density. Indeed this is so, provided we define the (regularized) ADM energy density of 
the background as 9 

EV 3 = lim {E r V 3 } = lim (-— ^" / V^t 2 K V3 dv 3 (r) ) , (6.23) 

where the 3-boundary v 3 (r) is the spacelike foliation of dAi$ and 2 K v:i is its extrinsic 
curvature. Explicit evaluation of (6.23) yields 

3 



E r = 



8ttG 5 



4A+B j^i 



(6.24) 



which is indeed the second term of (6.22). 

Having derived the general asymptotic expansions for the black hole and the PW 
geometry in (6.7) and (6.8), we are now ready to evaluate AF: 



Fbh — Fpw — — (Sbh — Spw) T + (Ebh — Epw) 



1 



1 



,4A 



dy_ 

Or 



horizon, BH 



— 3 lim 

r— >oo 



e 4A+B A , 



- e 4A+B A' 



BH 



PW. 



8ttG, 



-e 4a 5 - 3 A™ f , 



(6.25) 



where in the last line we have used (3.20). The evaluation of the limit in (6.25) is 
rather simple. We choose a direct matching condition of the BH and PW boundaries, 
parameterized by x [see (6.8) and (6.9)] and x [see (6.7)] respectively: 



x = So x . 



(6.26) 



Additionally we have to set 



B 



0. 



(6.27) 



PW 



Matching the boundary values of the scalars p and x fc> r the black hole and reference 
geometries yields 

Pii = I k 2 5l x = k 5 . (6.28) 
Furthermore, matching the asymptotic volumes of the BH and PW phases determines 

1 



So 



(6.29) 



3 As usual, the reference background has to be subtracted before the r — > oo limit is taken. 
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The final result is 



Apw = e 4? (-§ + i V 2 Pio -^V'ln^e V ^ ) , (6.30) 
where we have introduced 

77 = ifee" € . (6.31) 
The difference of free energies thus has the form 



p 4a / 

Fbh ~ F PW = -— U + 3 e" 4Q A^ 

;r ' V:i / -'") (5 + 3 , Aft,' 



(6.32) 



, 27r 

This can be further simplified by using the integral of motion, (3.19). Indeed, evaluating 
the constant in (3.19) in the IR and the UV, and equating them, we find 

ln/3 + ln2 + 4£ = 4a + m<5. (6.33) 

Thus we can rewrite (6.32) as 

F BH - F PW = --g^ ^ (ft + V % - 1 ^ In (± e, 2 )) . (6.34) 

Notice that from (6.8) the residual reparametrization invariance x — > Ax can be ab- 
sorbed by the following transformation on the quantities {£, Pio, Pn, Xo> Xio, 

Pio A 2 pi + A 2 pn In A , /On — > A 2 pn , 

4 (6.35) 
Xo^Axo, Xio -»■ >?Xio + -A 2 XolnA, 



3 



P -> A 4 /3 . 



This leaves (6.34) invariant. 

From the gauge theory arguments, we expect that the free energy of the PW phase 
scales as N 1 . On the other hand, the iV-scaling in (6.34) suggests that in the large 
TV-limit, F PW is essentially zero compared to F BH . Thus, in the high temperature 
phase, we identify the Af = 2* gauge theory Helmholtz free energy density F at large 
't Hooft coupling with (6.34) 

— 771 j-i — rp jrenomn 

= -r BH — fPW = 1 I E 

vr 2 iV 2 1 A , ( 1 ,\\ (6-36) 



4(2 



28 



Also, from (5.18) and (6.33), the N = 2* gauge theory entropy density S is 

9-9 - ^ e 3 «- ^ ( 2(3 Y\* (6 37) 

s = Sbh -8{^) ~W^{t) 6 ■ (6 - 37) 

Finally, the gauge theory energy density E is that of the (renormalized) ADM energy 
density [from (6.25)] 

E = E BH - E PW = T I™ + TS BH . (6.38) 

The gauge theory temperature T is identified with that of the horizon in the BH phase, 
(5.16). 

6.2 The high temperature thermodynamics of the J\f = 2* 

Given the analytical expression for the high temperature expansion, (5.5)-(5.14), it 
is straightforward to determine the leading correction to the non-extremal D3 brane 
thermodynamics due to the PW mass flow. As we will note, the sign of this correction 
suggests the possibility of the "deconfinement — > finite temperature Coulomb phase" 
phase transition. 

We begin by matching the {a,ai,a 2 } parameters in (5.5) with {£, x, p} of the 
asymptotic expansion, (6.8). Recall that the y coordinate in (5.10) is just e B . Thus 
the x coordinate of (6.8) and y are related according to 

y~ l-2x A . (6.39) 

This corresponds to setting (3 = 2 in (6.8). Now, to linear order in a,, by matching the 
scalars p, x m (5.10) and (6.8), we find 

41n2 8 

Pio — «i, Pu — «i, 

, 2[r(|)] 4 ( 6 - 4 °) 

Xo = TwsTW "2 , Xio " 



12 z ' ^ ±u n 2 



[r(i)F 

Notice that xio is independent of a 2 . Furthermore, matching to the asymptotic PW 
solution, (6.28) and (6.29), determines 



vr.o^ [r(i)] 2 



Notice that, to leading order, a,\ oc a\. This is consistent with the gauge theory 
expectation that, to leading order in the M = 2 hypermultiplet mass, it is enough to 
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turn on only the mass for the fermionic components. Thus the consistency of the high 
temperature expansion (conditional to the asymptotic Af = 2 supersymmetry) requires 
setting a\ to zero. 

Now, matching A, we find 

= a-§lii2 + a|| 2 -f, a = a + a%£ 2 - (6.42) 

where [using (5.13)] 



6 = A 



2/-1- 



with X2 given by (5.10). There is a nontrivial check on the computation. With (6.42) 
and (5.15), we find from (6.33) 

(In2) + ln2 + 4(d-|ln2 + a 2 2 | 2 ) = 4 (a + a 2 , f 2 ) + a 2 (I + | 72 ) , (6.44) 

or 

4(6-6)=| + |72 = 4^. (6.45) 

Given expression (6.43), we have numerically verified that (6.45) is indeed correct. 

Using (5.16), (6.36) and (6.37), we can now express the gauge theory thermody- 
namic quantities in the high temperature regime in terms of = ke~^ <C 1}: 



1 je (, , r(3/4) 



T = — — eM 1 + 



4 



2V27T \ 8vr 2 



r/ 



S=^e^l-I^, 2 ) (6.46) 
25/% 1 8tt 2 77 ' 



/V 2 

F = E-ST = -— -e 4 « 
32vr 2 



Recalling (6.4), and inverting the T <-> £ relation above 
we finally obtain the thermodynamic quantities 



^^(.-^^ (6.48, 
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For the thermodynamic process (6.48) we find 10 

TdS^dE. (6.49) 
In the next subsection we discuss a possible resolution of this puzzle, (6.49). 

6.3 The chemical potential of the J\f = 2* flow? 

We suggest here 11 that the apparent violation of the first law of thermodynamics for 
the leading in m/T C 1 correction to the high temperature thermodynamics of the 
M = 2* gauge theory, (6.49), could be explained as due to the neglection of the in- 
duced chemical potential for the temperature deformed M = 2* flow. We stress that, 
while a certain chemical potential appears to resolve the paradox, we do not have an 
understanding of what exactly is its corresponding conjugate operator. Additionally, 
it is conceivable that a different subtraction procedure for the computation of the su- 
pergravity effective action and the ADM mass would resolve the problem with the first 
law of thermodynamics altogether [28]. Having said this, however, here we restrict our 
attention to the possibility of having an induced chemical potential for the temperature 
deformed PW flow. 

One of the basic statements of the gauge/string theory correspondence [2] is the 
identification of the type IIB string theory partition function with the M = 4 gauge 
theory partition function, where the boundary values of the string fields $ are the 
sources of the gauge theory operators 



Z string]^] 



= Z gauge [<S> ] = e -"W[*o] , (6.50) 

<3>(i — >oo)=<J>o 



where W gauge is the generating functional for the connected Green's function in the 
gauge theory 

W gauge [* ] = -\n(eS**«> *) . (6.51) 



gauge [ 

\ I gauge 

It is not known how to precisely define the string theory partition function. But, ig- 
noring all the stringy a' corrections 12 , and also all the string loop corrections (which 

10 We have confirmed the leading correction to the free energy in (6.48), and thus the violation of 
the first law of thermodynamics, numerically. For details see section 6.4. 

11 We would like to thank Chris Herzog, David Lowe and Andrei Starinets for very useful discussions. 

12 Neglecting a' corrections implies that the string fields $ must actually be type IIB supergravity 
modes. 
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basically amounts to taking the N — > oo limit with the large but finite 't Hooft cou- 
pling), it is reasonable to assume that Z string is dominated by its saddle point 13 - 
the extremum of the (Euclidean) supergravity action Ie with the prescribed boundary 
values of the sources $o : 



In Z s tring[^} 



(6.52) 

*(r^oo)=<I>o 



extremum Ie[<&] 

>oo)=$oy 

We restrict to the Af = 4 gauge theory deformations which are irrelevant in the 
UV (the finite temperature Af = 2* flow discussed in previous sections is precisely of 
this type). This implies that the asymptotic geometry that extremizes Ie is necessarily 
AdS 5 x S 5 

ds 2 10 ~ e~ 2r/L dxl + dr 2 + L 2 dQ 2 5 . (6.53) 

r— >oo 

Generically, the supergravity mode $ that extremizes Ie behaves as 



$ e (A-4)r/L +JFo e -Ar/L > ( g_ 54) 



where A is the mass dimension of the gauge theory operator (9*, and Tq should be 
interpreted as its vacuum expectation value: 

(Off|0*|O H > =^ , (6.55) 

where |0h) is a vacuum state of the deformed Af = 4 Hamiltonian H 

H = H N=A + $ o $ . (6.56) 

As was emphasized in [19], in a theory with a unique (or at least isolated) vacuum, 
the dynamics should determine the vev (6.55) once the Hamiltonian (6.56) is specified. 
Generically we expect an isolated vacuum whenever the Af — 4 supersymmetry is 
completely broken. This will always be the case whenever arbitrary deformations of the 
type (6.56) are supplemented by the finite temperature deformation. Consider now such 
a deformation, namely gauge theory with Hamiltonian (6.56) at finite temperature. 
From the gauge theory perspective we can definite two different partition functions: a 
canonical partition function 14 

Z„[T] = e-^ = TTe^ H , (6.57) 



13 The subtleties of multiple saddle points will not arise in the present situation, namely the high 
temperature phase of the string theory dual to Af = 2* gauge theory 
14 We assume that the gauge theory volume is constant. 
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where the Tr is taken over the eigenstates of the full Hamiltonian, or the grand canonical 
partition function 



^gauge 



[T, na] = e -^ n[T '^ ] = Tr e -7 H +i**Q* . (6.58) 



There the trace is taken over the eigenstates of both the Hamiltonian H and the 
(conserved) "charge" operator Qq>, conjugate to the chemical potential 

When we neglect the fluctuations in Q$, we obtain 

n[T, (i< s> ]=-T In E gauge [T, n*] 

(6.60) 

~f[t]-^(q*)v 3 , 

where V 3 is the spatial volume of the gauge theory coming from the integration over 
the zero momentum modes. If we identify the string theory partition function (6.52) 
with the canonical partition function of the gauge theory (6.57), or equivalently 

F[T] = Tin /e±f d3x *°°\ ~ T extremum I E [<S>] 

\ I gauge 

in the case of the finite temperature M = 2* gauge/string duality we will face the 
breakdown of the first law of thermodynamics, (6.49). Rather, we propose that one 
should identify the string theory partition function (6.52) with the grand canonical 
partition function of the gauge theory (6.60), or equivalently, 

n[T,n*} =T\n/e±f d3x *°°\ = Tin / e -£ J* 3 * 

\ ' gauge \ I gauge 

~F[T}-n*(Q*)V 3 (6.62) 



(6.61) 

<E>(i — >oo)=<J> 



~ T extremum Ie[$\ 



= TI E [%], 

<E>(i — >oo)=$o 



where the last equivalence defines Je[$o]- Notice that the one-point correlation function 
can be computed by differentiating with respect to /i$ the correspondence (6.62): 

From (6.62) and (6.63) we find 



t(^J d 3 x(Q^ = -V 3 {Qz) . (6.63) 



F[T] = n[T,^]-^ dn[ ^^ ] . (6.64) 
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Finally, the first law of thermodynamics (for dV^—O) takes the form 

dF = -S dT + fi* d((Q*)) , (6.65) 

where the independent variables are //$ and T. 

In the rest of this subsection we demonstrate that with interpretation (6.62), there 
is no conflict with the high temperature thermodynamics of the M = 2* flow. First 
of all, notice that though the M = 2* flow necessarily has moduli, its finite tempera- 
ture deformation should not. Thus we expect that turning on the mass term for the 
fermions (supergravity dual to the five- dimensional scalar x) should uniquely fix their 
condensate. In the language of the asymptotic behavior of the scalar x m (6-8), this 
implies that specifying xo should uniquely determine the coefficient of its normalizable 
mode oc XoXio- This is precisely what we find in (6.40). If in our case we take $ = m, 
we would obtain \x = —m/T. Following (6.62) 15 , 

n[T,/i] = TI E [m\ = + £ m 2 T 2 = + £ ^ , (6.66) 

where in the second equality we have substituted the Helmholtz free energy from (6.48), 
which by computation equals T I E [m]. Additionally, to avoid cluttering the formulas 

we set 

2 



^^r(3/4) 4 . (6.67) 



7T 



Notice that the entropy is (6.48) 

vr 2 iV 2 T 3 K 



S = - ■- /i 2 T 3 . (6.68) 



From (6.63) and (6.64), we find 

(Q,) = -^ 4 , 

8 8 ^ 

It is easy to see that the first law of thermodynamics, (6.65), is now satisfied: 

._(^_^ r s) dr + /ld (_^) (6.70) 
= -S dT + fid((Q^)) . 



15 



As before we talk about densities of the thermodynamic quantities. 
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Thus we have shown that if we assume that the finite temperature deformation of 
the PW flow has an induced chemical potential p = —m/T, the interpretation of the 
supergravity computation in terms of the grand canonical ensemble appears to resolve 
the puzzle with the first law, (6.49). What is not clear, however, is what is exactly 
the charge operator conjugate to p. Though it appears that the expectation value 
of Q^, (6.69), is related to the gaugino condensate, Q M cannot be the fermion mass 
operator; the latter does not commute with the gauge theory Hamiltonian and thus 
cannot be conserved. 

6.4 The phase transition 

Independent of the high temperature expansion, the general expression for the gener- 
alized free energy density difference between the BH and the PW phases is given by 
(6.34): 

5Q[T, p] = n BH - n PW = -^1 e 4 « (p + V 2 p 10 - 1 r? 4 In Q e ^ ) . (6.71) 

Following the discussion of the previous subsection, we have reinterpreted the Helmholtz 
free energy F as the generalized free energy Q [T, p] . This expression is valid for arbi- 
trary temperature T and chemical potential p = —m/T, which in turn are implicitly 
related to the parameters of the supergravity solution, £, f3, rj and p w , which show up 
on the right hand side of (6.71). 

To proceed beyond the high temperature expansion, we may examine the behavior 
of SQ numerically. To do so, we extract the appropriate coefficients governing the 
behavior of Sfl by matching the UV behavior of the numerical solution with (6.8) and 
the IR behavior with (3.20). In particular, we first work in the UV and fix the x- 
coordinate of (6.8) through the functional dependence of the scalar B (recalling that 
(3 may be scaled away). Then, after matching the coefficients of the leading nontrivial 
asymptotics {f>n, Xo} with the asymptotic PW geometry according to (6.28) and (6.29), 
we may unambiguously extract the subleading terms {pio,Xio}- We finally obtain £ 
through the relation (6.33), where a and S are determined from the behavior of A and 
B at the horizon. 

Of course, {£,PiO)Xio} are functions of the data at the horizon {a,Po ? Xo}? (3.20). 
Actually p an d Xo cannot be independent since the coefficients of the leading UV 
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Figure 3: Numerical determination of the generalized free energy 5fl (6.71) with 
the leading behavior 5Qq = — 7r 2 iV 2 T 4 /8 subtracted. The numerical values are 
given by the data points, while the leading high-temperature correction, 5Q\ = 
r(3/4) 4 iV 2 m 2 T 2 /87r 2 , i s indicated by the dashed line. 



asymptotics of p and x-> {PiijXo}, must satisfy (6.28 

Pn 2 



2 • ( 6 - 72 ) 



which is just the statement of asymptotic M = 2 supersymmetry. This results in a 
reduction to two parameters, T and fi (or equivalently T and m). Finally, since any 
scale in the M = 2* theory may be related to m, we note that, for the numerical work, 
we only need to examine a one parameter set of solutions. 

For a thermodynamic process at a fixed volume, the physical phase is realized from 
the minimization of the generalized free energy 

n[T, = min {n BH , n PW } . (6.73) 

Thus the signature of a phase transition would be the vanishing of Sfl at a certain 
critical temperature T c 



Sfi [T, n = -m/T) 



= 0. (6.74) 

T=T C 
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In the high temperature phase, m/T <C 1, we have found [compare with (6.66)] 



that is, SQ < 0. On the other hand, we have argued that in the low temperature phase, 
m/T >> 1, we would instead expect 5Q > 0; see Fig. 1. 

The result of the numerical work is shown in Fig. 3. While the numerics appear 
to be in good agreement 16 with our analytical prediction (6.75), we have been unable 
to confirm the phase transition. One possibility is that the critical temperature of the 
conjectured phase transition is at T c = m g, where g is a small number. This would 
make numerical study of the transition rather challenging, as the ultra-low temperature 
supergravity flows (see Fig. 2) approach the supersymmetric (singular) PW flow, and 
are plagued by numerical instabilities. Another possibility is that in the large N 
limit this phase transition is actually at T c = 0. This issue clearly deserves further 
investigation. 
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